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Abstract 

The energy dissipated by viriahzation shocks during hierarchical structure formation of the Galaxy can 
exceed that injected by concomitant supernova (SN) explosions. Cosmic rays (CRs) accelerated by such 
shocks may therefore dominate over SNe in the production of ^Li through a + a fusion without co-producing 
Be and B. This process can give a more natural account of the observed ^Li abundance in metal-poor stars 
compared to standard SN CR scenarios. Future searches for correlations between the ®Li abundance and 
the kinematic properties of halo stars may constitute an important probe of how the Galaxy and its halo 
formed. Furthermore, ®Li may offer interesting clues to some fundamental but currently unresolved issues 
in cosmology and structure formation on sub-galactic scales. 
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1 Introduction 

Apart from ^Li, the bulk of the Hght elements Li, Be and 
B are believed to arise from nonthermal nuclear reactions 
induced by cosmic rays (CRs) (Meneguzzi , Audouze & 
Reeves 1971). In the last decade, extensive observations of 
LiBeB in population II, metal-poor halo stars (MPHS) have 
turned up new and unexpected results, spurring controversy 
as to what type of CR sources and production mechanisms 
were operating in the halo of the early, forming Galaxy (e.g. 
Vangioni-Flam, Casse & Audouze 2000). To date, most 
models of hght element evolution in the early Galaxy have 
focused on strong shocks driven by supernovae (SNe) as the 
principal sources of CRs. A general consensus is that Be 
and B in MPHS mainly originate from the "inverse" spal- 
lation process, whereby CR CNO particles are transformed 
in flight into LiBeB by impinging on ISM H or He atoms 
(Duncan, Lambert & Lemke 1992). This can be realized 
if a sizable fraction of the CRs responsible for spallation 
comprise fresh, CNO-rich SN ejecta (e.g. Vangioni-Flam 
et al. 2000, Ramaty et al. 2000, Suzuki, Yoshii & Kajino 
1999, Suzuki & Yoshii 2001, hereafter SY), as opposed to 
CRs injected from the average ISM (e.g. Fields & Olive 
1999). 

The origin of ''Li in MPHS (Hobbs 2000 and references 
therein) is more mysterious, as current models involving SN 
CRs face some difficulties. A peculiar aspect of Li is that in 
addition to spallation, the fusion process of CR a particles 
with ambient He atoms can be effective, and should ac- 
tually dominate Li production at low metallicities. (Note 
that while both ^Li and ®Li are synthesized in compara- 
ble amounts, the CR-produced ^Li component is generally 
overwhelmed by the "Spite plateau" from primordial nucle- 
osynthesis in the metallicity range under consideration; e.g. 
Ryan et al. 2001.) If the CR energy spectrum is taken to be 
a standard power-law distribution in momentum (§3), one 
requires a CR injection efficiency much higher than nor- 
mally inferred to reproduce the ®Li observations, whether 
the CR composition is metal-enriched or not (Ramaty et 
al. 2000; SY, see their flg.2). This raises the question of 
whether there may have been other sources for ®Li. 



Suzuki & Inoue (2002; hereafter SI) proposed a new and 
more natural ''Li production scenario based on a previously 
unconsidered CR source: CRs accelerated at virialization 
shocks (VSs), i.e. shocks driven by the gravitational infall 
and merging of sub- Galactic gas clumps during the hier- 
archical build-up of structure in the early Galaxy. Such 
shocks are inevitable consequences in the currently stan- 
dard theory of hierarchical structure formation (see how- 
ever §5). We further discuss a number of testable predic- 
tions and important implications expected in our scenario 
for understanding the formation of our Galaxy. The po- 
tential relevance to some fundamental issues in cosmology 
and structure formation on sub-galactic scales is also men- 
tioned. 



2 Cosmic Ray Sources in the Early Galaxy 

SNe are known to release Esn ~ 10^^ erg of kinetic en- 
ergy in each explosion, driving strong shocks into the am- 
bient medium. These SN shocks are favorable sites for effi- 
cient CR acceleration through the first order Fermi mecha- 
nism. A plausible value for the injection efficiency fsN, i.e. 
the fraction of the SN kinetic energy imparted to CRs, is 
10 — 20%, deduced from comparison of the SN rate and the 
energy content of CRs currently observed in the Galactic 
disk. The global energetics of SNe for the early Galactic 
halo can be estimated from the cosmic star formation his- 
tory or the total amount of heavy elements ejected by halo 
SNe (see SI for details of the estimation). Both methods 
give esN ~ £sNl-imp/AIg ~ 0.15 keV per particle as the 
average specific energy input from SNe. 

SNe are not the only sources of mechanical energy (and 
hence CRs through shock acceleration) that may have been 
active in the early Galaxy. In the framework of the cur- 
rently successful picture of hierarchical structure formation 
in the universe, large scale objects are formed through the 
merging and virialization of smaller subsystems, driven by 
gravitational forces acting on the dark matter. For each 
merging hierarchy, shocks should inevitably arise in the as- 
sociated baryonic gas component, whereby the kinetic en- 
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ergy of infall is dissipated and the gas heated to the virial 
temperature of the merged halo (e.g. White & Rees 1978). 
It is quite plausible that such virialization shocks (VSs) also 
accelerate CRs (e.g. Miniati et al. 2001). 

A guide to how structure formation may have proceeded 
in our Galaxy, particularly for the Galactic halo, may be 
offered by the recent numerical simulations of Galaxy for- 
mation by Bekki & Chiba (2001, hereafter BC; see also 
Samland & Gerhard 2003, Abadi et al. 2003). It is ex- 
pected that sub-Galactic structure eventually merge into 
a single entity in the central region at redshift z ~ 2, 
whereby the majority of the infall kinetic energy is viri- 
alized. The total energy dissipated at this main VS can be 
evaluated from the virial temperature of the merged sys- 
tem. The virial temperature T„ for a halo of total mass 
Mt virializing at redshift z ~ 2 is /cbT^ = fj,mpGMt/2rv — 
0.26keV (Mt/3 x IQ^'^Mof^^ , where r„ is the virial radius, 
and we have assumed a cosmology with Q,m. — 0.3, JIa = 0.7 
and h = 0.7 (see SI). The specific energy dissipated at the 
VS is thus evs ~ 3fc_gr„/2 ~ 0.4 keV per particle, higher 
than the above estimate for SNe by a factor of ~ 2.6, if 
we adopt A/t = 3 x 10^^ (Sakamoto, Chiba & Beers 
2003). The CR contribution from VSs compared to SNe 
should similarly be higher, as the CR injection efficiency 
^vs should not be too different from that for SNe (Miniati 
et al. 2001). CRs accelerated by VSs should therefore be 
at least as important as SN CRs, and may well dominate 
at early epochs. 

Another, more speculative but potentially interesting 
possibility is outflows powered by massive black hole(s), 
which may have been active in the early Galaxy. However, 
there are numerous ambiguities with such a picture, and 
here we will concentrate on CRs from VSs and SNe (see SI 
for more discussion). 



3 Model 

We employ assumptions and parameters deemed most plau- 
sible for the SN CRs. The CR injection efflciency is taken to 
be 5sN = 0.15 for each SN of kinetic energy _Esn = 10^^ erg. 
A standard, single power-law distribution in particle mo- 
mentum is adopted for the CR spectrum (see SI). The in- 
jection spectral index is chosen to be 7sn = 2.1, appropriate 
for strong SN shocks, and consistent with the source spec- 
trum inferred for present-day CRs. As assumed in SY, the 
composition of SN CRs is a mixture of SN ejecta contain- 
ing freshly synthesized CNO and Fe and the ambient ISM 
swept up by the SN blastwave. 

CRs from VSs are markedly different from SN CRs in 
a number of important ways. First and foremost, VSs do 
not synthesize fresh CNO nor Fe, so that the composition 
of these CRs is completely ascertained by the pre-existing 
ISM. When the ISM is metal-poor, these shocks induce very 
httle Be or B production through inverse spallation, and 
are only efficient at spawning Li via a — a fusion. Second, 
VSs are not necessarily strong ones, particularly for major 
mergers of systems with comparable masses (Miniati et al. 
2001), provided that the pre-shock gas has not cooled sig- 
nificantly. Shock acceleration should then lead to injection 
indices 7vs steeper than the strong shock limit value of 2, 
which works in favor of Li production (§4). Major merger 



shocks may possess Mach numbers as low as ~ 2 — 3, cor- 
responding to 7vs — 2.5 — 3.3; 7vs = 3 is generally chosen 
below. As with SN CRs, we take the spectral shape to 
be a momentum power-law distribution and the injection 
efficiency to be ^vs = 0.15. 

A further distinction from SNe is that the VS CR flux 
should not entail any direct dependence on the metallicity, 
which is in fact an obstacle to predictive modeling. While 
the hierarchical growth of structure with respect to cosmic 
time may be evaluated in concrete ways using e.g. extended 
Press-Schechter formalisms, relating this to [Fe/H] requires 
additional knowledge of how the combination of star for- 
mation, SN nucleosynthesis and chemical evolution in our 
Galaxy proceeded with cosmic time. This involves large 
uncertainties, and is not specified in our chemical evolution 
model. As a first step, we choose to describe the time evolu- 
tion of VS CRs in a simple, parameterized way, assuming a 
'step function' behavior: VS CRs begin to be injected from 
a certain time fys, maintains a constant fiux for a duration 
TVS, and then returns to zero. We take the injection du- 
ration TVS to be roughly the dynamical time of the major 
merger, ~ 3 x 10* yr. The injection is also assumed to be 
uniform, since the effect of the main VS should be global 
throughout the gas under consideration, unlike SNe. The 
injected VS CR flux integrated over rvs is normalized to 
the above values for ^vs and evs- The true evolutionary 
behavior of VS activity relative to metallicity should actu- 
ally be probed through future observations of ^Li at low 
[Fe/H] (§4, §5). 

After injection by either SNe or VSs, the spectral flux 
Fi{E,t) for each CR element i evolves with time during 
subsequent interstellar propagation. This is obtained from 
time-dependent solutions of the CR transport equation for a 
leaky box propagation model (SY). Using the transported 
spectra, we calculate the CR production of LiBeB in the 
ISM including all three types of reactions: forward spalla- 
tion of ISM CNO by CR protons and a's, inverse spallation 
of CR CNOs by ISM H and He, and the fusion of CR a's 
with ISM He. For more details, consult SY. We do not 
consider here the potentially complicating effects of stellar 
depletion, which are highly uncertain at the moment. 



4 Results 

For our calculations, we have selected the following sets 
of parameters for tvs, ^vs and 7vs, respectively, labeled 
models I - IV: I (0.22, 0.1, 3), H (0.22, 0.1, 2), HI (0.32, 
0.1, 3), and IV (0.1, 0.5, 3), where t\rs and rvs are in units 
of Gyr. These were chosen to provide results exemplary of 
light element production by VS CRs, in contradistinction to 
that by SN CRs. The evolution of ® Li and Be vs. metallicity 
calculated for each model until the end of halo chemical 
evolution ([Fe/H]~ —1.5) is shown in Fig0 along with the 
current observational data in MPHS. 

We discuss some salient points regarding these results. 
First, it is conflrmed that with our flducial parameters, pro- 
duction by SN CRs alone (dashed) works very well for the 
observed Be (and B, not shown), yet falls short of the ob- 
served ®Li. In contrast, with reasonable values for evs, ^vs 
and 7vs , production by VS CRs is capable of explaining the 
current ®Li data quite adequately. This mainly owes to two 
facts: 1) VSs are more energetic than (or at least as ener- 
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Figure 1: Model results of ^Li/R (thick) and Bc/H 
(thin) vs. [Fe/H], for SN CRs only (dashed curves) and 
SN plus VS CRs (solid curves), each label correspond- 
ing to the parameter set described in the text. Also 
plotted are current observational data for ^Li (thick 
crosses from Smith et al. 1998 and Nissen et al.2000; 
squares with errors from Asplund et al. 2001; circle 
with upper limit from Aoki et al. 2003) and Be (thin 
crosses from Boesgaard et al. 1999). 

gctic as) SN shocks, as estimated in §2, and 2) VS CRs can 
generate ®Li at early epochs independently of the metallic- 
ity. Regardless of the early evolutionary behavior, identical 
®Li abundances are attained at the end of the halo phase 
for a given 7vs (I, HI, & IV), since this is determined by the 
time-integrated CR flux, for which we had assumed a fixed 
value. Compared to a flat spectral index of tvs = 2 (II), 
a steeper one of tvs = 3, more appropriate for low Mach 
number VS shocks (§3), results in a larger ®Li yield, by 
about factor of 3. This is because with a constant total CR 
energy, a steeper index implies a larger CR flux in the sub- 
relativistic energy range E ^ 10 - 100 MeV, where the ®Li 
production cross section peaks (see fig.l in SY). Note that 
a conservatively lower specific energy for VSs, evs ~ 0.15 
keV/particle (i.e. comparable to SNe), can still bo consis- 
tent with the available ®Li data provided that 7vs — 3. 

As already stressed, 2) is a consequence of ®Li synthe- 
sis being dominated by a — a fusion, and VSs not creating 
any new Fo. Depending on the onset time and duration of 
the VS, the ''Li abundance may potentially reach largo val- 
ues quickly at very low metallicity, which can be followed 
by a plateau or a very slow rise. On the other hand, SNe 
unavoidably give forth to freshly synthesized Fe, so a cor- 
relation with Fo/H must arise; in fact ''Li/H vs. Fc/H for 
SN CRs can never be much flatter than linear (e.g. Ra- 
maty et al. 2000). Moreover, since VSs do not eject fresh 
CNO, they produce very little Be or B through spallation; 
only a minuscule contribution can appear as the ISM be- 
comes metal-enriched. This may allow an extremely large 



^Li/Bc ratio at low [Fe/H]. Conversely, we see that SN CRs 
must play an indispensable role in generating the Be and B 
observed in MPHS. 

Independent of the particular evolutionary parameters, 
the following abundance trends are characteristic of VS CR 
production and should serve as distinguishing properties of 
the scenario for future observations. Going from high to low 
metallicity: a plateau or a very slow decrease in log®Li/H 
vs. [Fc/H], followed by a steeper decline in some range of 
[Fc/H] corresponding to the main epoch of VS; a steady 
increase in ''Li/Be, possibly up to values exceeding ~ 100, 
also followed by a downturn. These traits are very dis- 
tinctive and not expected in SN CR models, for which the 
slope of log''Li/H - [Fe/H] must be ~ 1 or greater, and the 
®Li/Be ratio constant at sufficiently low [Fe/H]. Distinc- 
tion from any production processes in the early universe 
(e.g. Jedamzik 2000) should also be straightforward, as 
they predict a true plateau down to the lowest [Fe/H], in 
contrast to an eventual decrease for VS CR models. Fur- 
ther, unique diagnostic features are discussed in §5. 



5 ''Li as Fossil Record of Dissipative Processes 
during Galaxy Formation 

A truly unique and intriguing cispect of the VS CR picture 
is that ®Li in MPHS can be interpreted and utilized as a 
fossil record of dissipative gas dynamical processes in the 
early Galaxy. Of particular interest arc various correlations 
expected between the ®Li abundance and the kinematic 
properties of the stars. On the one hand, ®Li arises as a 
consequence of gaseous dissipation through gravitationally- 
drivcn shocks, and survives to this day as signatures of the 
dynamical history of hierarchical structure formation in the 
early Galaxy. On the other, the kinematic characteristics of 
stars presently observed should reflect the past dynamical 
state of their parent gas systems, because once stars form, 
they become collisionloss and have long timcscalcs for phase 
space mixing (e.g. Chiba & Beers 2000). Interesting rela- 
tionships may then exist among the two observables. 

For example, recent intensive studies of the structure 
and kinematics of MPHS in our Galaxy based on Hippar- 
cos data (Chiba & Beers 2000) have elucidated the de- 
tailed characteristics of our Galaxy's halo, such as its two- 
component nature: an inner halo which is flattened and 
rotating, and an outer halo which is spherical and non- 
rotating. This dichotomy has been suggested to result from 
diff(>r(:>ii(;c;s in the physical processes responsible for their 
formation, dissipative gas dynamics being crucial for the 
former, and dissipationless stellar dynamics determining 
the latter. The numerical simulations of Galaxy forma- 
tion by BC support this conjecture: the outer halo forms 
through dissipationless merging of small sub-Galactic clumps 
that have already turned into stars (c.f. Searle & Zinn, 
1978), whereas the inner halo mainly forms through dissi- 
pative merging and accretion of larger clumps that are still 
gas rich (c.f. Eggen, Lynden-Bell & Sandage, 1962). In our 
scenario, ^Li production is a direct outcome of the principal 
gas dissipation mechanism of gravitational shock heating. 
If the above inferences on the formation of halo structure 
arc correct, ^Li should bo systematically more abundant 
in stars belonging to the inner halo compared to those of 
the outer halo, a clearly testable prediction. An important 
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prospect is that ®Li may provide a quantitative measure 
of the efTectiveness of gas dynamical processes during for- 
mation of halo structure, rather than just the qualitative 
deductions allowed by kinematic studies. 

Another possibility regards the main epoch of VS with 
respect to metallicity (§3). As already mentioned, the re- 
lation between ®Li/H and Fe/H should mirror the time 
evolution of dissipative energy release through VSs, but is 
complicated by being convolved with the uncertain ingre- 
dients of star formation and chemical evolution. The Fe 
abundance can be a bad tracer of time, especially at low 
[Fe/H] where the effects of dispersion in SN yields can be ex- 
tremely large (SY). Stellar kinematics information may of- 
fer a handle on this problem. The observed relation between 
[Fe/H] and < >, the mean azimuthal rotation velocity of 
MPHS, seems to manifest a distinctive kink around [Fe/H] 
~ -2 (Chiba & Beers 2000). Through chemo-dynamical 
modeling of the early Galaxy, BC have proposed that this 
kink may correspond to the epoch of the major merger (§3). 
If this was true, a simple expectation in the context of the 
VS CR model is that ®Li/H should be just rising near this 
value of [Fe/H], which is the range occupied by the cur- 
rently ^Li-detected stars; also expected are a steep decline 
at lower [Fe/H], as well as a plateau or slow rise at higher 
[Fe/H] (i.e. close to model HI in fig0. However, any infer- 
ences related to [Fe/H] are always subject to the chemical 
evolution ambiguities. A more reliable and quantitative an- 
swer may be achieved by looking for correlations between 
^Li/H and < > without recourse to [Fe/H], as ^Li is 
a direct and pure indicator of dynamical evolution in the 
early Galaxy. 

Thus the VS CR model for ®Li bears important im- 
plications for understanding how our Galaxy formed. If 
the above mentioned trends are indeed observed, it would 
not only confirm the virialization shock origin of ®Li, but 
may potentially point to new studies of "^Li Galactic ar- 
chaeology" , whereby extensive observations of ®Li in MPHS 
coupled with detailed chemo-dynamical models can be ex- 
ploited as a robust and clear-cut probe of dissipative dy- 
namical processes that were essential for the formation of 
the Galaxy. 

Furthermore, ®Li in Galactic MPHSs may potentially 
offer interesting clues to a number of outstanding current 
problems in cosmology and structure formation theory, all 
involving physics on sub-galactic scales. 1) The global im- 
portance of dynamical feedback by SNe has been a long- 
standing uncertainty in galaxy formation theory (White & 
Rees 1978, Abadi et al. 2003). Further observations of ''Li 
vs. Fe or possibly Be and B at low metallicity and com- 
parison with detailed theoretical models may constrain this 
crucial unknown. 2) The efficiency and location of VSs as 
commonly assumed on sub-galactic scales has recently been 
brought into question (Katz et al. 2002, Birnboim & Dekel 
2003), with radical implications for how stars and galax- 
ies form and how the galaxy luminosity function is shaped 
(Binney 2003). The early evolution of ''Li may provide 
a direct probe of this presently speculative but important 
suggestion. 3) Comparison of theoretical simulations with 
observations of dwarf galaxy cores and of satellite galaxies 
of our Galaxy and within the Local Group indicate that 
standard cold dark matter (CDM) produces much more 
substructure than is actually seen. This "CDM crisis" may 
point to non-standard dark matter properties, such as warm 



dark matter, self-interacting dark matter or even more ex- 
otic proposals (e.g. Ostriker & Steinhardt 2003, Madau & 
Kuhlen 2003). On the other hand, the apparent discrep- 
ancy may be the result of strong feedback by SNe or a UV 
background. These different possibilities should result in 
differences in the early evolution of ''Li, potentially consti- 
tuting a unique probe. 4) From a combined analysis of mea- 
surements of cosmic microwave background anisotropies by 
WMAP and of the power spectrum on galactic and sub- 
galactic scales, it has been suggested that the spectrum of 
primordial density fluctuations deviates from a standard, 
single power-law (Spergel et al. 2003). Although less dras- 
tic than non-standard dark matter, this would also modify 
the growth of structure on subgalactic scales (Madau & 
Kuhlen 2003), and the consequent ®Li production at early 
epochs. More detailed investigations of these intriguing 
prospects are certainly necessary, but in principle, ®Li in our 
Galactic halo may shed light on these issues of paramount 
importance for cosmology and the physics of the early Uni- 
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